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CALCULATION  OF  GROUND  THAWING  ALLOWING  FOR  WATER  SEEPAGE 


S trol tel 'stvo  v  Rayonakh  Voatochnoy  Slblri  G.  M.  Fel'dman 

1  Kraynego  Severe  (Construction  in  Regions 
of  Eastern  Siberia  and  Far  North),  Collec. 

12,  Krasnoyarsk,  1969,  pp  111.144 


As  Is  known,  the  nature  of  the  heat  exchange  of  the  filtration  flow 
with  frocen  soil  depends  significantly  on  the  filtering  properties  of  this 
soil  In  a  frosen  state  [l"|.  We  should  differentiate  two  possible  cases* 
the  soil  pores  are  filled  with  Ice  or  not  (open  pores). 

Open  Pores  (Fig.  1).  On  the  soil  surface,  there  is  maintained  a  layer 
of  water  H  with  temperature  Tj.  The  y.axts  Is  directed  downward;  the  origin 

of  coordinates  Is  on  the  soil  surface.  Up  to  the  time  7  ,  the  frosen  soil 
has  thawed  to  depth  .  The  filtration  factor  and  the  porosity  of  thawed 
soil  (0  v  j  <  V)  Is  signified  by  Kj;  ,  while  that  of  frosen  soil  Is 
1  nd  1  cated  by  (  >_  ^  y<vy)-k2;m2» 


At  the  time  moment  {  ,  the  water  has  seeped  to  depth  Jr  .  Let  us  assume 

that  during  the  filtration  of  water.  In  the  frosen  soils,  the  plugging  of 
pores  will  not  occur.  Let  us  further  assume  that  at  each  time  moment  7, 
the  temperature  distribution  In  the  thawed  ground  corresponds  to  the  estab¬ 
lished  state.  Moreover,  the  velocity  v  of  water's  velocity  will  depend  only 
on  time. 


The 
System  of 


determination  of  function  £,(  f)  reduces  to  a  solution  of  the 

equations  for  heat  conductivity  and  filtration* 


a£t_vKs0 

adyl  V  dl/ 


v  <  y< 


dy, _ **  _ _ H»  y« 
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(2) 


2 


where  /l  •  the  heat  conductance  coefficient  of  thawed  soil,  kcal/m.hr°Cj 
T 

C  ‘ “  heat  capacity  and  volumetric  weight  of  wateri  v  -  velocity  of  water 
miveaent,  m/hr|  and  t  -  soil  temperature,  °C, 

A  solution  to  this  system  of  equations  for  yQl  Vjt;  as  a  function  of 

1  does  not  appear  possible.  However,  with  an  adequate  degree  of  approxi¬ 
mation,  we  can  recotsnend  the  following  procedure  for  determining  £  -  t.  ( 't  ). 


The  solution  to  Eq,  (1)  has  the  forms 

_  t.-t.  '-** 


t-f iafrfc-iK 


(5) 


v-z-\ 

The  solution  of  Eq,  (2)  at  t  -  const  has  the  form  {^2 ”)s 

1/  _  1/  H*  Vd 


li  f  %\()  y**»4 M  ^  \f  -  k  — LLtJ^ - 


(6) 


where  A  -  (Kj/Kj  -  1)  \  /H. 

If  we  assume  that  Kj  -  K2 ;  m^  -  ,  Eq.  (6)  will  acquire  a  simpler 

form  [2")s 


(7) 


Substituting  Eq.  (5)  Into  the  Stefan  condition  (4),  we  derive 

a(^-D  - 

After  Integration  of  this  expression,  we  derives 

‘^(£’5-0  <»> 

Assigning  the  time  value  ^j,  based  on  Eq.  (7),  we  find  the  corresponding 
values  for  y^  and  v^.  Then  substituting  the  ^  and  Vj-values  into  Eq.  (8), 
we  find  the  value. 


This  method  of 


calculating  twill  be  sufficiently  accurate  becauses 


r  ~\J  f  1  dT  nf 


For  the  following  time  value  ?2»  the  values  of  y_0 , 
(6),  substituting  Into  this  formula  ?- 
-  -?2  and  v  -  Vj  Into  (8)  or  (9),  we  find  the  value 
At  |vjf/it|  >  3,  Eq.  (8)  can  be  wrlttens 


V2  will  be  found  with 
Further,  substituting 
for  and  so  forth 


Eq. 

f- 
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(9) 


*  cx  .  Ar(t|-tj)rv.^ 

f'V  csvvo  ' 

Closed  Pores.  In  a  case  of  closed  pores,  the  necessity  originates  of 
calculating  the  conduct  of  heat  exchange  between  the  filtration  flow  and  the 
subjacent  permafrost  soils  (Fig.  2). 


Me  place  the  origin  of  coordinates  on  the  soil  surface.  The  positive 
y-axls  Is  directed  downward.  On  the  soil  surface,  temperature  t  -  t2«  At 
time  moment  1  ,  the  depth  of  thawing  fc,  j  at  this  depth,  t  -  0.  The  direction 
of  filtration  flow  Is  along  the  X-axls.  The  original  water  temperature  at 
X  -  0  equals  t  -  tj. 


A  simitar  problem  has  been  solved  by  V.  G.  Gol 'dtmap  [lj  for  determin¬ 
ing  the  effective  heat  conductivity  factor  (In  the  direction  of  the  Y- 

axls)  and  also  the  thawing  depth  of  frozen  groundi 


Cldv*  *  or  • 
t(v,o)rt,  ; 
I  ^  n  : 
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u  y  iv 
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(10) 


As  Is  evident  from  the  diagram  (10),  In  a  calculation  of  the  conductive  heat 
exchange,  we  Incorrectly  assigned  the  boundary  condition  at  the  upper  surface 
of  filtration  flow..  it/^y|^_0  -  0.  In  reality,  on  the  upper  surface  of 
filtration  flow,  we  should  not  assume  the  absence  of  a  heat  flow  but  a  quite 
definite  temperature  value  t  -  t2  equalling  the  soil  temperature  under  natural 


.4. 


conditions  at  the  occurrence  depth  of  the  filtration  flow's  upper  surface 
(In  Fig.  2,  y  -  O).  In  addition,  we  did  not  take  Into  account  the  conductive 
transport  of  heat  alonR  the  filtration  flow.  A  more  precise  mathematical 
formulation  of  the  stipulated  problem  has  the  forms 


3\t_  _  V  at 

T>  y”  o 
y  -  0.  t  -tf; 

VaY.t=0- 


ot 


..  Jf-Oi  o<VOf;t»t, 


V 


(ID 


where  v  -  velocity  of  water  movement,  m/hri 

Cl  3  »  m2  /hr 


/lef  m  effective  heat  conductivity  factor  In  direction  of  y-axls; 
c  and  -  heat  capacity  and  volumetric  weight  of  water,  respectively. 


Figure  3. 

A  solution  to  the  problem  has  the  forms 


-5- 


vJ  /-  ‘ ”75  U'  /  Pl  vy»4 

Let  us  determine  the  heat  flow  at  the  lower  boundary  y  -  i 


T  •  o.| 


The  value  of  the  sums 


where 


r  n.l  ,  l  1 

\  [(-l)  »  l]c.x?l-  *  V  V’-  j 


<*  v  v 


Is  represented  In  Fig.  3.  As  Is  obvious  from  the  figure,  the  value  for  the 
sue  can  be  represented  In  the  formi 

B-if:  jK-»$£H«p[-  “*> 


where  the  constants  A,  B  and  C  are  selected  on  the  basis  of  the  curve  (Fig. 


3)  In  the  known  range  of  variations  In  criterion  v  - 


If  the  water  Is 


filtered  by  the  complete  section  through  a  soli  layer  which  has  thawed,  the 
velocity  of  water's  movement  equals i 

V  •  J7jT  ^  (15) 

where  Q  -  output  of  water,  kg/hour;  L  -  width  of  flow,  m;  and  tp  ■  soil 
temperature  on  surface,  °C.  If  the  upper  surface  of  filtration  flow  occurs 
at  depth  t  from  the  day  surface.  It  follows  thati 

■f.  "til'1?  ) _  U; - — — -  (16) 

I  V  Lr(?_£)  (16) 

The  effective  heat  conductivity  factor  according  to  V.  G.  Gol'dtman  consti¬ 
tutes  a  function  of  velocltyi 

Aj-e*  V&  (17) 

where  Am  heat  conductivity  factor  of  water-saturated  thawed  coarse-grained 
soli;  0  -  coefficient  of  heat  dispersion.  According  to  V.  G.  Gol'dtman,  for 
the  Magadan  region,  D  -  7,  Let  us  analyse  the  case  of  a  complete  section. 
Substituting  (14),  (15)  and  (17)  Into  (13),  we  derives 


Let  us  formulate  the  Stefan  equation  at  the  boundary  of  phase  transitional 
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where  Cm  heat  of  phase  transitions  kcal/k2;  and  W  •  volumetric  Iciness  of 
frosen  soil,  k2/m3.  Dividing  the  variables,  we  derives 


V 

GW 


1 


(20) 


Substituting  the  "a**  value  into  (20),  we  gets 

f,  t  dr 


v_ 

csw 


(20) 


For  the  low  x-values,  the  function  ^5  ■  (>(x^  2  )  can  be  determined  by  a  numer¬ 
ical  integration  of  Eq,  (20).  For  the  high  values,  Eq.  (20)  will  assume  the 
forms 

WC  ^  ?■  dir 


G  W 


1  4-  I* »— 
*  Lf* 


(21) 


After  Integration  and  transformations,  we  obtain  the  following  expres¬ 
sion  for  estimating  the  thawing  depth  at  high  X-valuess 


t>a 


AtjT.T1 _ 

crv/  "  2  /ILf 


iM-fJ- 


ua 


(22) 


The  system  (10)  can  be  utilised  for  estimating  the  thawing  rate  of  two 
lateral  sides  of  a  vertical  tallk  crevice  In  which  a  filtration  flow  will 
move,  or  the  base  and  top  of  a  water-bearing  layer  bounded  by  frosen  layers. 
The  mathematical  formulation  of  such  a  problem  has  the  forms 

}  ,  \  _  r  3t_  _  q  "J 

+A«F5»  7*  l 

•  J(>0;  t-O;  J  (23) 

?; 


j<  rO]  yet1?;  trt« 


Moreover,  in  the  given  case  is  half  of  the  thawed  layer,  M.  We  match  the 
origin  of  coordinates  with  the  middle  of  the  thawed  soil  layer.  The  solu¬ 
tion  of  problem  (23)  has  the  forms 


t-tfjt-0  ro5(£^-)^F{^2Ap)1"  xf i  ^ 


(24) 


where  -  (2n-l)  *1/2 ;  n-  1,  2,  3.... 
boundary  of  phase  transltlonsi 


We  determine  the  heat  flow  at  the 


T  n»«  ^ 


(25) 


Now  let  us  solve  a  similar  problem  under  the  condition  that  we  disregard  the 
conductive  heat  transfer  in  the  direction  of  the  main  filtration  flow  (X-axis) 
as  compared  with  the  convective  flow  in  this  same  direction. 


Then 


lit  _  v_iL„ot 

fT1  o  d* 
j(  i>0;  y*  t  Y ;  t  =0j 

^0;V<^,  t*  t, 

where  2  !;■  depth  of  thawed  soli  layer,  m.  The  solution  to  the  problem  has 
the  forms  t 


On 


Let  us  determine  the  heat  flow  at  the  boundary  of  the  phase  transitions 


(27) 


Let  us  formulate  the  Stefan  condition  for  estimating  the  thawing  rate  of 
frozen  grounds 


0A.  t.t'  [  <i[2n-0‘!'jll  df 

>7vT~J  ’ G  w  dr 


The  velocity  of  water  movement i 


V  =  7T 


a 


°lt^  ' 

t  of  el 


(23) 


(29) 


By  analogy  with  (17),  the  coefficient  of  effective  heat  conductivity  X  f  has 
the  forms  e 


Coefficient  "a" 


n  .2XLIft5S 
2  L-,r  T  e 


Substituting  (29)-(3l)  Into  (28),  after  transf or mat  1  on ,  we  derives 

ZXS.  Jr _  r-d* 


kii  .( - - - 


(2q.i)M(2X  Lr7»T>QtJ 
Ma^c  y1 


*] 


(30) 


(31) 


(32) 


where  2  V,  ■  Initial  soli  depth,  m.  The  estimation  of  the  depth  of  thawed 

t  O 

8 ol 1  2  V.  can  be  conducted  by  a  numerical  Integration  of  Eq,  (32). 

On  a  calculation  of  aclcular  hydro-thawl ng.  The  mathematical  formula¬ 
tion  of  the  problem  has  the  forms 

1  at  .  _  vyr  .0;  , 

A*<t 

OeO  ;  i  ;  t  *t,; 


/ 


(33) 
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The  following  notations  have  been  adopted!  h  -  depth  of  needle's  submergence, 
m|  tj  -  Initial  water  temperature;  v  -  velocity  of  water  flow,  m/hr;  a  - 

ef/c  I  c  and  Y  -  heat  capacity  and  weight  of  water  by  volumo;  ?-  radius, 
m;  radius  of  halo  of  soil's  thawing  around  needle,  m;  lf  -  effective  co¬ 

efficient  of  heat  conductivity,  kcal/m.hr,  °C;  and  A  -  heat  conductivity  co¬ 
efficient  of  water-saturated  thawed  ground,  kcal/m.hr,  °C. 


The  solution  of  problem  (33)  has  the  form! 

where  fu  ^  -  roots  of  equation  JQ  (h)  -  0. 

-2/.05;  >i»8,rj4,  >v-11.8;  >/- M,9j 

A  A“2,*2i  A-24*85;  A"5*32- 


(34) 
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I,  I?  *<lnJ  jfrL  the  Bessel  functions  of  the  flr3t  type,  of  zeroth  and  first 
orders,  respectively.  We  determine  velocity  v  at  assigned  output  Q  kg/hri 


-  <)/>.  i,2, 


We  then  finds 


ex z. 


05) 


(36) 


We  find  the  value  for  the  heat  flow  to  the  boundary  of  phase  transitions 


(with  consideration  of  Eos,  (35)  and  (36)). 


DO. 


~ ( '^  nr  r »A 


u) 


(37) 


Let  us  examine  a  case  wnen  cne  value  of  conductive  heat  flow  along  the  X-axls 
can  be  disregarded.  Then  Eqs.  (34)  and  (37)  will  acquire  a  simpler  forms 


(38) 


(39) 


The  value  of  the  sum  ^Texrjl  *r depending  on  parameter  Tl^ax/Q  has  been 

presented  In  Fig.  4.  As  Is  apparent  from  the  figure.  In  a  fixed  range  of 
variation  In  parameter  \  j ax/Q,  the  sum's  value  can  be  presented  In  the  forms 


(40) 


i-i 


where  A,  B  and  C  •  parameters  of  the  curve  (see  Fig.  4).  Substituting  (36) 
and  (40)  Into  (39),  we  obtains 


(41) 


Let  us  write  the  Stefan  equation  for  estimating  the  thawing  halo  of  frozen 
ground  around  a  needle  V,(Xj,/l)s 


^  ( kCi 


(42) 


Separating  the  variables,  we  get  the  following  equations 

f 

2t,T 


G  w 


(43) 


where  7M  p  -  needle's  radius,  m.  If  we  assume  a  -  acp»  !#*•  adopt  the  average 
value  of  the  conventional  temperature  conductivity  factor  for  the  entire  period 
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of  hydro- thawing,  Eq.  (43)  wlil  assune  the  forai 


&*j[oc£jr 

a  + 


(44) 


Integrating  Eq.  (44),  we  get  the  following  expression  for  determining  the 
thawing  halo  of  the  soil  around  a  needle i 


iffZlrbQ 


(45) 
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